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^ ■ Abstract 

' We discuss the generalized von Neumann (Tsallis) entropy and the gener- 

ahzed Fisher information (GFI) in nonextensive quantum systems, by using 
Q ■ the interpolation approximation (lA) which has been shown to yield good 

1^ i results for the quantal distributions within 0{q — 1) and in high- and low- 

temperature limits, q being the entropic index [H. Hasegawa, Phys. Rev. E 
80 (2009) 011126]. Three types of GFIs which have been proposed so far in 
the nonextensive statistics, are discussed from the viewpoint of their metric 
properties and the Cramer-Rao theorem. Numerical calculations of the q- 
' and temperature-dependent Tsallis entropy and GFIs are performed for the 

electron band model and the Debye phonon model. 
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1 Introduction 



The Boltzmann-Gibbs-Shannon-von Neumann entropy and the Fisher information 
play central roles as information measures in classical and quantum statistics (for 
review see [I], relevant references therein). The Boltzmann-Gibbs-Shannon-von Neu- 
mann entropy represents a global measure of ignorance, while the Fisher information 
expresses a loal measure of a positive amount of information pj. The Fishser infor- 
mation has the two properties: (i) it expresses the metric tensor for the neighboring 
points in the Riemannian space spanned by probability distributions or density ma- 
trices, and (ii) it provides the lower bound of unbiased estimation errors in the 
Cramer- Rao theorem (CRT). The Fisher information has been employed for a study 
of effecientcy of information transmission and its decodings. 

In the last decade, much progress has been made in the nonextensive statistics 
initiated by Tsallis [2], who proposed the generalized entropy (called the Tsallis 
entropy) defined by 



k 



B 



q-1 



for classical case, 



for quantum case. 



(1) 
(2) 



Here q denotes the entropic index, ks the Boltzmann constant, Pq{x) the probability 
distribution and Pg the density matrix. The Tsallis entropy is a one-parameter 
generalization of the Boltzmann-Gibbs-Shannon-von Neumann entropy, to which it 
reduces in the limit of g 1.0. The Tsallis entropy is non-additive in a sense that 
for p{A UB) = p{A) p{B) or p{A U B) = p{A) ® p{A), we obtain 



SgiAUB) = SgiA) + SgiB) 
^ Sg{A) + Sg{B) 



k 



B 



-SM)SgiB), 
for q ^ I. 



(3) 
(4) 



The Tsallis entropy is super-extesnsive (sub-extensive) for g < 1 (g > 1), and 
q — 1 expresses the degree of the nonextensivity. The nonextensive statistics has 
been widely applied to various subjects in physics, chemistry, information science, 
biology and economics [3j. 

At the moment, three types of generalized Fisher informations (GFIs) have been 
proposed in the nonextensive statistics [1]-[10], 



Gf) 



q Tr 



Tr 



Tr 




(5) 
(6) 

(7) 



q y— Pq/Tr p^) expresses the escort density matrix and 6'„ a parameter 



where 

specifying pq-. the classical case of Eqs. 



is obtainable if we read pq — > Pq{x) 
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and Tr ^ j dx [see Eqs. (l85!) - (l87j) ]. In the limit of g — 1.0, the three GFIs given 
by Eqs. ([5]), and ([7j) reduce to the conventional expression, 



The GFI in Eq. (E]) [H [5] which is derived with the use of the generalized KuUback- 
Leibler divergence [HI [121 [13] , expresses the metric tensor of the first properties (i) 
of the Fisher information discussed above. The GFI in Eq. ([5]), however, is not 
applicable to the CRT for variance and covariance of physical quantities averaged 
over the escort density matrix. The GFI given by Eq. ([6]) [6, 7j preserves the second 
properties (ii) of the lower bound in the CRT, although it does not have the metric 
properties (i). The GFI given by Eq. (E]) [8], [9l has both the properties (i) and 
(ii). Details of the three GFIs and a comparison among them will be discussed in 
this paper (Sees. Ill and V). 

The Tsallis entropy and generalized Fisher information (GFI) in nonextensive 
c/asszca/ systems have been considerably studied in [H [6], El [3 [10] [H]- [19]. In recent 
years, the generalized von Neumann {quantum Tsallis) entropy has been investigated 
in the bipartite spin-1/2 [20], Id spin-1/2 [21j, 2d bosonic systems [21] and Hubbard 
dimers [22]. The purpose of this paper is to discuss the Tsallis entropy as well as the 
GFI in nonextensive quantum systems. Physical quantities in nonextensive quan- 
tum systems are evaluated by the trace over the escort density matrix [23], which 
may be formally expressed in exact integral representations [2H [25] . Their actual 
evaluations are, however, tedious and difficult because they involve self-consistent 
calculations of expectation values of the energy and number of particles. This is the 
case in a calculation of the generalized Bose-Einstein and Fermi-Dirac distributions 
(referred to g-BED and g-FDD, respectively) in nonextensive quantum statistics. 
Quite recently it has been pointed out that this difficulty may be overcome when 
we adopt the interpolation approximation (lA) which yields good results within the 
0{q — 1) and in the high- and low-temperature limits [26j. Indeed the g-BED and 
g-FDD in the lA are given in simple analytic expressions, which have been suc- 
cessfully applied to nonextensive quantum systems such as black-body radiation, 
Bose-Einstein condensation, BCS superconductivity and metallic ferromagnetism 
^27] . The three GFIs which are formally expressed in exact integral representations, 
may be much simplified when we adopt the lA [26J. 

The paper is organized as follows. In Sec. II, we briefly explain the exact and 
interpolation approaches to nonextensive quantum statistics [26]. In Sec. Ill, the 
Tsallis entropy and the three GFIs in the lA are discussed. Numerical calculations 
are reported in Sec. IV for the electron band model and for the Debye phonon model. 
Discussion and conclusion are presented in Sec. V, where a comparison among the 
three GFIs is made. The CRT for the three GFIs in [H [Qj [TOl [7] is discussed also for 
the g-Gaussian distribution in the nonextensive classical statistics. 
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2 Exact and interpolation approaches 
2.1 Exact approach 

We first obtain the optimum density matrix p^, applying the maximum entropy 
method (MEM) with the optimum Langange multipher (OLM) [23] to the general- 
ized von Neumann (Tsallis) entropy given by Eq. ([2]) under the constraints given 

by 

Trp, = 1, (9) 

Tr{plN} = TrplN,, 

Tr{plH} = TrplE,, 

where Ng and Eg denote expectation values of the number operator (A^) and the 
Hamiltonian (H), respectively. The OLM-MEM leads to the density matrix given 
by [23] 

1 , , 



with 



Xg = Trw, (11) 
Ng = -^Tr{wm}, (12) 

Eg = -^Tr{wm}, (13) 



w 



[1 - (1 - q)(3{H -pN-Eg + pNg)Y/^'-'>\ (14) 



where f3 stands for the inverse temperature and p the chemical potential (Fermi 
level). The escort density matrix Pg is expressed by 

p" 1 

Pg = ^ = ^W^, (15) 

Cg Aq 



where we have employed the relation, 



c, 



Trpl = Xl~^. (16) 



By using the formulae for the gamma function T{z) given by 



CO 



= —— / n^-^e"^" du for s > 0, (17) 

r s 



x' = — r(s + 1) / (-t)"""^e-^*rft for3fts>0, 

^TT Jc 



we may express Eqs. ( ITT|) - (fT3|) as integrals along the real axis for q > 1.0 and in the 
complex plane for g < 1.0 as follows [211 [251 126] : 

/~ G (u; ^, du for g > 1, 

^'^^ ^/^(^;t^'(T^)^iH)^^ forg<l, ^^^^ 

c 



i: /o~ G FT- 7^) >l(w)iVi(M)c/M for g > 1, 
N„= { > ^ < (20) 

^ 2^jH[t;^^,j^)U-t)N^i-t)dt forg<l, 

c 

i^/o~^(^;FT' (^)^iH^iM^^ forg>l, 



where 



Y^{u) = e"(^^-'^^')Si(M), (22) 
E,{u) = ^re-"(^-^*) = e^S^'■^[^^^""^^'="''^ (23) 
Niiu) = 5^/i(efc,w), (24) 

k 

Ei{u) = J2 fii^k,u)ek, (25) 

k 

/^(^'") = e"(-^)Tl ' ^^^^ 

G(n;a,6) = -^n'^-ie"^", (27) 
r(a) 

H{t;a,b) = T{a + l)b-'' i-ty-^e-^'. (28) 

Upper and lower signs in Eqs. (l23l) and (l26l) are applied to boson and fermion, 
respectively, and C denotes the Hankel path in the complex plane [211 125] . Note 
that Nq and Eg are obtained by Eqs. (fT9|) - (!2T]) in a self-consistent way. Such self- 
consistent calculations have been made for the electron band model and the Debye 
phonon model in [26] . 

2.2 Interpolation approach 

Self-consistent calculations including Nq and Eg are rather tedious and difficult. In 
order to overcome this difficulty, we have proposed the lA [26], assuming that 



Xq Xq 



e^(^^-t^^^^) Ei{u) = 1, (29) 
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in Eqs. (!20l) and (!2T!) . Then they are expressed by 



u 



1 



' (g-l)/3 



27r 



1 



-g' (i-g)/3 



A''i(u) for q > I, 
Ni{-t) dt for g < 1, 



(30) 



Equations fl20l) and fl2T|) are alternatively expressed by 



(31) 



(32) 
(33) 



where the quantal distribution /q(ejt,/3) is given by 



c 



i^lH(t;j^^,j^,)fi{e,,-t)dt forg<l. 



(34) 



With the use of Eq. flMl) . the analytic expression of the g-BED in the lA is given 
by [26] 



/.(e'/5) = ^[^^"^"^'^1' for0<g<3, 



?i=0 



where x = /3(e — /i) and expresses the g-exponential function defined by 



= [1 + (1 - q)x\ 



1/(1-9) 



(35) 



(36) 



with [?/]+ = y Q{y) and 6(?/) is the Heaviside function. 

Similarly, the analytic expression of the g-FDD in the lA is given by ^26] 



for e > /i, 
1 for e = /i, 

1.0 - Fq{\x\) for e < /i, 



(37) 



with 



{n+l)x-\q 



for < g < 3, 



(38) 



n=0 



where x = j3{e — fi). // (e, (3) given by Eqs. (!35l) and (j37l) reduces to /i(e, /?) in the 



limit of g ^ 1.0 where — > e^. 



Sg = kB^—f = kB\ng{X,), (39) 



3 Information Entropies 

3.1 Tsallis entropy 

By using Eqs. ([2]) and f|T6l) . we obtain 

Q 

where Xg is given by Eq. (fT9l) . or alternatively by 

f rr G f -r^) e"S^(^^~'^)-^''(^^'")±Sfc''^[i±^i(^'")] rfw, for g > 1 , 
= <^ V ' 9-1' (3-i)/3y ' ^ ' (40) 

2-K JC^^ l-g' (l-g)/3 ) ^ iui y ^ i , 

Inq(x) being the g-logarithmic function, 

ln,(a;) = (41) 
In the limit of g — > 1.0 where Ing(x) — > ln(x), Eq. (HOl) reduces to 

which yields the well-known expression of the quantum Boltzmann-Gibbs entropy 
given by 

5i = A;Bln(Xi) = fcB/5(^(e-/i)/i(efe,/3)±5^[l±/i(efc,/?)]), 

\ k k J 

= -fcB ($^/iln/iT$^(l±/i)ln(l±/i)) . (42) 

\ k k ) 

3.2 Generalized Fisher information matrix 
3.2.1 ge^e^ 

The distance between two operators p and a in the Riemann space is defined by 

D,{p\\^) = Kgip\\^)+Kgia\\p), (43) 
with the generalized KuUback-Leibler divergence i^g(/3||(j) [TT| I 



K,{m = Y^iTrp-Trip'^a'-'^}), (44) 

which is in conformity with the Tsallis entropy given by Eq. ([2]). The distance 
between the neighboring operators of Pg({6'„}) (= Pg) and Pq{{On + 69n}) (= p'g) is 
expressed by 

DMp'g) = 5^ 5^^^™, (45) 



where the GFI of 5'6»„e„, is given by [5J 



d In pq\ f dlnpg 



the bracket {Q)q denoting the expectation value of an operator Q over p^, 

(Q), = Tr{pqQ}. (47) 



The GFI given by Eq. (H6|) has a clear geometrical meaning expressing the metric 
between the adjacent density matrices in the Riemannian space spanned by the 
OLM density matrices of [pq] [5J. Equations (jSj) and ( H6|) are quantum extensions 
of the counterparts in the nonextensive classical statistics 
For {61,62) = i/3,/3p), Eqs. and (06]) yield 



9ii 

922 
912 



921 = - Tr {w'^-\H - Eq){N - iV,)}. 



(48) 
(49) 
(50) 



3.2.2 Gf) n 

Naudts proposed the quantum GFI given by [7] 



Gf) 



-idpq 



d6r. 



-idpq 



d6r^ 



where the bracket [Q]q denotes the expectation value given by 



Q 



Tr {PqQ}. 



(51) 



(52) 



When we consider the expectation value of [A\q for an operator of A and 6n = 6m = 6, 
the generalized CRT is shown to be expressed by [7J 

fdA'\2 



U - [A] 



> 



(53) 



where 



A' 



{A)q. (54) 

Note that A' is an average over pq. The CRT given by Eq. ( I53ll may be derived 
as follows [?]. Taking the derivative of A' with respect to 6 and using the relation: 
dpq/d6 = 0, we obtain 

dA' „ ^dpq 
'dd 



d6 



Tr{^-^{A-[A]q)}, 



Tr{Pq[Pq'"^-^]{A-[A]q)}. 



(55) 
(56) 



Employing the Cauchy-Schwartz inequality, we obtain 



dA' 



< Tr{P,{p-'^] }Tr{P,{A~[A],r}, 



(57) 



which leads to Eq. ( l53l) . 

It is straightforward to extend the method mentioned above to the case of 
(Ai,y42) = {H,N) and {9i,92) = {[3,(3^). The generalized CRT is expressed by 



V > C"^ C = D 

where the covariance matrix V is given by 
/ 



(58) 



V 



V 



{H-E,){N-N,) 



{H-E,){N-N,) 



J 



Calculations using Eqs. f|T5|) and fl5T]) yield elements of the Fisher information 
matrix G given by 



(H - E,f 



Gil = 
G22 = 
G12 = G21 = 

and those of C expressed by 

Gil -- 
G22 - 
G12 - 
G21 - 

A simple calculation leads to 



11; 



22, 



iH-E,)iN-N, 



12 



-V. 



21, 



(59) 
(60) 
(61) 



9(/?/i) 
d{(3fi) 
d(3 



H{H - Eg 
N{N - Ng) 
H{N - Ng 
N{H - E, 



-I Q 



-I 9 



D = V, 

which implies that the CRT given by Eq. (158|1 is satisfied with an equal sign. 



(62) 
(63) 
(64) 
(65) 

(66) 
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3-2.3 ge„e„. 

The GFI given by 



d6m. 



(67) 



provides the lower bound of unbiased estimates in the CRT as shown in the following. 
For the expectation value of A = [A]q, we obtain 



dA ^ , / J r n \ -1 

_ . Tr{^iA-[A],)}, 



(68) 
(69) 



By using the Cauchy-Schwartz inequality, we then obtain the CRT given by 



U - [A] 



(dA\2 
\2i ^ \ de ) 



> 



(70) 



With the use of Eq. ( |T5l) and ( 1671) . the generalized CRT for the case of {Ai, A2) 
{H,N) and (6*1,^2) = is expressed by 



V > c"^ C = D, 

where elements of g are given by 



9n 

922 



X, 



Tr {w"^-'iH-E,Y}, 



912 = ^2i = -(£) Tr{^3'^-2(^-E,)(iV-iV,)}, 



(71) 

(72) 
(73) 
(74) 



and those of C are given by 

Cu - 
C22 - 
C12 - 
C21 - 







d(3 




ON, - 




dm ^ 








dm ^ 








d/3 





^Tr{w^'^-^H{H-E,)}, 



(75) 
(76) 
(77) 
(78) 



It is noted note that the GFI given by Eq. (l67l) is nothing but a quantum extension 
of that proposed in nonextensive classical statistics [8l 191 flO]. 
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Table 1: A comparison among three GFIs for On = Om = [q ^ 1.0) 



GFI 


Metric 


CRT 


Refs. 


gge = qTr{p-\dpJd9f} 


Eq. (05]) 


NA 


mm 


Gee = Tr{p-\dpjdef} 


NA 


Eq. dSSl) 




~gee = Tr{p-\dPjdef} 


Eq. (Ci) 


Eq. (CO]) 





CRT: Cramer-Rao theorem, NA: not apphcable 



We may show that ge^Om given by 
neighboring escort density matrices of 
as given by 

Di{P,\\P'^ = 
where Di{p\\cr) stands for the distance 



Cq. (I^7|) denotes the metric tensor for the 
\i{en}) {= P,) and P,({^„ + 56^}) {= P'^) 

Yl ~9s^o^ ^(^n59m, (79) 

nm 

;iven by Eqs. and dB]) with q = 1.0. 



3.2.4 A comparison among g, G and g 

In the preceding subsections, we have discussed the properties of g, G and g. A 
comparison among the three GFIs is made in Table 1. We note that g expresses the 
metric in the Riemann space spanned by density matrices (pq), but G does not. In 
contrast, G provides the lower bound for the CRT, to which g is not applicable. On 
the contrary, g denotes the distance between the escort density matrices (Pg) and it 
also satisfies the CRT. It is shown that G provides a better bound for the inequality 
in the CRT than g: V = D > D. 

When we adopt the exact transformation with the use of formulae for the gamma 
function given by Eqs. (fTTI) and (fT8|) . equations for the GFIs may be expressed as 
integrals along real axis for g > 1.0 or in the complex plane for q < 1.0. Expressions 
of the GFIs for q > 1.0 with the lA [Eq. fl29l) ] are summarized in the Appendix. It 
is now possible for us to numerically calculate the Tsallis entropy Sg and the GFI 
matrices as functions of q and temperature. We will report such numerical calcu- 
lations for the electron band model and the Debye phonon model in the following 
section. 

4 Numerical calculations 
4.1 Electron band model 

We employ a band model for electrons with a uniform density of state given by [26] 

pie) = il/2W)eiW-\e\), (80) 
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where W denotes a half of the total band width. We adopt = 0.5, for which 
/i = 0.0 independently of the temperature because of the adopted symmetric density 
of states given by Eq. (IHOj) . 

Figure [T] shows the temperature dependence of Eg for various q calculated self- 
consistently with the use of Eqs. (IT^ - fl^ . With increasing q from unity, Eg at 
higher temperatures is decreased, although that at lower temperatures is increased 
(see also the inset of Fig. 1 of Ref. |10|). We have calculated the average energy also 
by using Eq. fl3Tl) in the lA. The ratio of E^^^^/Eg (= A) is plotted in the inset of Fig. 
[H For q = 1.5, for example, this ratio is changed from unity at low temperatures 
{kBT/W ~ 0.005) to about 0.991 at high temperatures {kBT/W ~ 1.0). For 1.0 < 
q < 1.5, an agreement between Eg and Eg^"^^ is much better. 

The temperature dependence of Sg for various q is shown in Fig. El With increas- 
ing q from unity, the temperature dependence of Sg at low temperatures becomes 
more significant but its saturated value at high temperatures become smaller. 

Sohd, chain and dotted curves in Fig. [3] show temperature dependence of Ga, gu 
and gu {i = 1,2), respectively, for g = 1.1: results for g = 1.0 are plotted by dashed 
curves. Note that Vu = Gu [Eqs. (l59l) and f l60l) ] which implies that the temperature 
dependence of the variance Vu is the same as that of Gu- All the GFIs show a similar 
temperature dependence. A closer inspection, however, shows that there are some 
differences between them. In particular, with increasing q from unity, gu and gu are 
increased while Gu is decreased. 

4.2 Debye phonon model 

We adopt the Debye model whose phonon density of states is given by [26] 




where A = 9Na/w^, Na denotes the number of atoms, u the phonon frequency and 
uj£) the Debye cutoff frequency. 

Figure m shows the temperature dependence of Eg for various q (with /i = 0.0), 
which is self-consistently calculated by Eqs. (fT9|) - (pTl) . It shows that Eg is larger for 
larger q. The inset of Fig. Hlshows the ratio of Eg^/ Eg (= A) where Eg^ is calculated 
within the lA [Eq.([3l])]. Although the ratio is not good at 0.02 < T/Td < 0.5, 
it becomes better at 1 < T/Td < 2, where To signifies the Debye temperature 
(fc^To = Itujd)- 

Figure O shows the temperature dependence of Sg. We note that for larger q, the 
temperature dependence at low temperatures becomes more steep and its saturated 
value at high temperatures becomes smaller. 

Solid, chain and dotted curves in Fig. [6] show the temperature dependence of 
Gu, gu and gu, respectively, for q = 1.1: results for g = 1.0 are shown by dashed 
curves. Note that the temperature dependence of the variance Vu is the same as that 
of Gu [Eq. (l59l) ]. All the GFIs show a similar temperature dependence, although 
with increasing q, magnitudes of Gu and gu becomes larger while that of gu becomes 
smaller. 




for < cij < ujo, 



(81) 
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Filled and open marks in Fig. [7|sliow the q dependence of D and D, respectively, 
for the Debye phonon model: note that D = V for the GFI in [7]. It is shown that 
the ratio of D/V, which is unity for q = 1.0, is gradually reduced with increasing q. 



5 Conclusion and discussion 



It is instructive to make a comparison among g, G and g for the g-Gaussian distri- 
bution in classical nonextensive statistics. For given mean (yUg) and variance (a 



the g-Gaussian distribution Pq{x) and its escort distribution Pq{x) are expressed by 

m 



Pqix) 



l-(l-g) 



vZa 



q) 



2z/(t2 



i/Ci-?) 



(82) 



(83) 



with 



q-1 ^ \ q-1 2 



27rcr„ 



for 1 < g < 3, 
for g = 1, 
for < g < 1, 



where i' = (3 — g)/2 and B{a,b) stands for the Beta function. The GFIs given by 
Eqs. ( H6l) . ( l5Ti) and ( 1671) in the classical case are expressed by 



9ij 
G 



q / pg{x) 



1 dpq{x) 
Pq{x) ddi 
1 dpq{x) 



9ij 



^'^"""Pqix) de, 

1 dPqjx) 

'^"^^^Pqix) 89, 



1 dpq{x) 
Pq{x) dOj 

1 dpg{x) 

Pq{x) 06, 
1 dPqjx) 

Pq{x) 89, 



dx. 



dx. 



dx. 



(85) 
(86) 
(87) 



Averages over Pq{x) and Pq{x) are expressed by {Q)q and [Q]q, respectively. 

For (Ai, A2) = {[x]q, [{6x)^]q) and (6'i, 62) = {fiq, a"^) where 5x = x — fiq, we obtain 

m 



\ 



(3-g) 

4<74 



g(9+i) 



(9+1) 
4(2g-l)<74 
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d[x]q d[x 



q 



d[{Sxf]q d[{Sxf]q 

The CRT is then expressed by 

' (3-g)(2q-l)ag Q 

V > c-g-c .5=1 ^(^-^) 







On the contrary, for the GFI of G [F, 7] with (A;, = ((x)^, ((^x)^)^) and 
h,02) = (yUg,ag), we obtain 

' 4 

G = I 

(3-9)2 







4{5-3g)a 



f _ I duq da^ 

d{{Sxf), d{{5x)^)q 





(l \ 







The CRT is expressed by 

V > C"^G-^C = D 



-I 







(5-3g) 



which shows that the CRT is satisfied with an equal sign. 

Sohd curves in Figs[8]^a) and[8]^b) show Du/a^ and 022/(7^, respectively, dotted 
curves expressing Vu/a"^ and V22/o"g- Chain curves in Figs [8](a) and [8](b) show 
Dii/a"^ and -D22/c"g, respectively. For a comparison, we show by the dashed curve in 
Fig. [8](b), l/(722<7q [10] ■ These figures clearly show that G and g preserve the CRT 
and that g is not applicable to the CRT because V22 < l/5'22 for ? < 1 in Fig- [Hl(b). 
The q dependence of D/V in Figs. [H](a) and[H](b) is similar to that for Debye phonon 
model shown in Fig. [71 

Adopting the factorization approximation (FA) to a calculation of the grand- 
canonical partition function, Biiyiikkilic, Demirhan and Giilec [28] derived the gen- 
erahzed quantal distribution given by 
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where the upper (lower) sign is apphed to g-BED (g-FDD). It has been shown that 
fq^i^k) may be alternatively derived by applying the variational condition to the 
entropy S^^ given by [29] 

= -E{[/f^(^'^')]^lnJ/f^(e.)]T[l±/f^(e.)]'^lnJl±/f^(e,)]}, (89) 

k 

with the constraints: 

k 

k 

Quite recently the g-BED and g-FDD given by Eq. flHHl) in the FA are criticized 
based on the exact approach [2Bj. It has been pointed out that the 0{q — l)-order 
contribution in the FA does not agree with that of the exact approach and that its 
g-FDD yields inappropriate results even qualitatively. This criticism is applied also 
to the expression for the FA entropy given by Eq. (l89l) . 

To summarize, we have discussed the generalized von Neumann (Tsallis) entropy 
and the GFI in nonextensive quantum systems, by using the lA Numerical 
calculations of the g- and temperature-dependent information entropies have been 
performed for the electron band model and the Debye phonon model. A comparison 
among the three GFIs (Table 1) has shown that for the CRT in the nonextensive 
statistics, we have to employ G [7J or g [Hl|9l|10] rather than g [5] which has a geomet- 
rical meaning derived from the generalized Kullback-Leibler divergence. Although 
our present discussion has been confined to g > 1.0 for nonextensive quantum sys- 
tems, it is necessary to extend our study to the case of g < 1.0, which is our future 
subject. 
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A The generalized Fisher information in the lA 



By using formulae for the gamma function given by Eqs. (ITTl) and (fT8|) and adopting 
the lA given by Eq. (129|) . we obtain following expressions for the GFIs. 
Elements of g in Eqs. P8l) - (!50l) for g > 1.0 are expressed by 



911 



G u 



g-1 '(g-l)/3 



E'i'\u)-2EqE,{u)+E' 



(Al) 
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922 



q G iu; 



912 = 921 = -q G [u 



2q-l 1 
q-l'{q-l)P 

2q-l 1 



N[^\u)-2NgNi{u) + N', 



with 



N['\u) = [N,iu)]' + J2fii^k,u)[l±Me,,u)l 

k 

Ef\u) = [i?iH]2 + 5^A(6,,«)[l±/i(6,,n)]e^ 

k 

d'^\u) = E^{u)N,{u) + Y,fli^k,u)[l±f^{ek,u)]ek, 



,2 
-ki 



(A2) 

(A3) 

(A4) 
(A5) 
(A6) 



where the upper (lower) sign is apphed to boson (fermion). Relevant results for 
g < 1.0 may be obtainable with a proper modification. 

Elements of matrix of G in Eqs. f l59|) - fl6T]) are expressed by 



Gil 
G22 
G12 



with 



^2 
HN 



G21 



G\^u- 
g(u- 



-Et 



N, 



1 ' 



HN 



1 



G u- 



q-V{q-l)P 

q 1 

q-r{q-l)(3 

q 1 

g-l'(g-l)/3 



Ef\u) du, 
D^iu) du. 



(A7) 
(A8) 

(A9) 

(AlO) 
(All) 
(A12) 



Elements of matrix of g in Eqs. ( 172l) -( 17^ are expressed by 

3g-2 1 



9ii 



922 



9l2 



q G [u; 



q G \u; 



921 = -q G { u] 



g-l'(g-l)/? 
3g- 2 1 

g-i'(g-i)/9 

3g-2 1 



[E^^\u) - 2EgEi{u) + E^g] du, 



(A13) 



[N^^\u) - 2NgNi{u) + N^] du. 



X 



g-l'(g-l)/?, 
[D^^\u) - EqNiiu) - N,Ei{u) + E,N,] du. 



(A14) 



(A15) 
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In the limit of g — 1.0, all the GFIs reduce to 

gfl2 = ^21 = 912 = — — ■ (A18) 
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Figure 1: (Color online) The temperature dependence of Eg of the electron band 
model for q = 1.0 (dashed curves), q = 1.1 (chain curves), q = 1.2 (dotted curves), 
g = 1.3 (solid curves) and q = 1.5 (double-chain curve), the inset showing the ratio 
of A = E'//Eg. 



Figure 2: (Color online) The temperature dependence of the Tsallis entropy Sq of 
the electron band model ior q — 1.0 (the dashed curve), q — 1.1 (the chain curve), 
q — 1.2 (the dotted curve), q — 1.3 (the sohd curve) and q — 1.5 (the double-chain 
curve) . 



Figure 3: (Color online) The temperature dependence of the component {i = 
1,2) of G (solid curves), g (chain curves) and g (dotted curves) for q = 1.1 of 
the electron band model, results for q = 1.0 being plotted by dotted curves for a 
comparison. 



Figure 4: (Color online) The temperature dependence of Eg of the Debye phonon 
model for q = 1.0 (dashed curves), q = 1.1 (chain curves), q = 1.2 (dotted curves), 
q = 1.3 (solid curves) and q = 1.5 (double-chain curves), the inset showing the ratio 
of A = E'//Eg. 



Figure 5: (Color online) The temperature dependence of the Tsallis entropy Sg of 
the Debye phonon model for q — 1.0 (the dashed curve), q — 1.1 (the chain curve), 
q — 1.2 (the dotted curve), g = 1.3 (the sohd curve) and q — 1.5 (the double-chain 
curve) . 



Figure 6: (Color online) The temperature dependence of the {i, i) component (f = 
1,2) of G (solid curves), g (chain curves) and g (dotted curves) for q = 1.1 of 
the Debye phonon model, results for q = 1.0 being plotted by dashed curves for a 
comparison. 



Figure 7: (Color online) The q dependence of Du (filled circles), D22 (filled squares), 
D12 (filled triangles), Dn (open circles), D22 (open squares) and D12 (open triangles) 
with T/Td — 1.0 for the Debye phonon model: note that Vij — Dij — 1, 2). 
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Figure 8: (Color online) The q dependence of (a) Vu/a'^ (the dotted curve), Dn/a"^ 
(the sohd curve), Du/a"^ (the chain curve) and {l/gua"^) (the dashed curve), and 

(b) that of 1^22/0-^ (the dotted curve), D22/(^q (the solid curve), D22/(^q (the chain 
curve) and l/g22(Tq (the dashed curve): note that Vu — Du — l/gu in (a) and 
V22 = -D22 in (b). 
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